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Abstract 

We calculate analytically quasi-normal modes of AdS Schwarzschild black holes includ- 
ing first-order corrections. We consider massive scalar, gravitational and electromagnetic 
perturbations. Our results are in good agreement with numerical calculations. In the case 
of electromagnetic perturbations, ours is the first calculation to provide an analytic expres- 
sion for quasi-normal frequencies, because the effective potential vanishes at zeroth order. 
We show that the first-order correction is logarithmic. 
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1 Introduction 



Recently there has been a lot of research studying quasi-normal modes of black holes in 
asymptotically AdS space-times [1]. Understanding these modes may give some insight into 
the AdS/CFT correspondence. 

Here we present a fairly comprehensive study of quasi-normal modes of AdS Schwarzschild 
black holes with a metric in d dimensions given by 

ds* = -f(r)df + *L + r*dnl_ 2 , f(r) = * + 1 - ^ . (1) 

and derive analytical expressions including first-order corrections. The results are in good 
agreement with results of numerical analysis. 

In the case of massive perturbations, which we discuss in section 2, we extend the ap- 
proach of [2] to include black holes of arbitrary size. We perform an expansion in 1/m, where 
m is the mass of the perturbation. The calculation involves a large amount of cancellations 
between various terms resulting in a sensible perturbative expansion of asymptotic expres- 
sions for quasi-normal frequencies. In section 3, we discuss gravitational perturbations and 
obtain the first-order corrections to the analytic expressions derived in [3, 4] (adapting the 
monodromy argument proposed in [5] and extended to first order in [6]). We find good 
agreement with numerical results [7]. In section 4, we extend the discussion to electromag- 
netic perturbations. In this case, the zeroth-order effective potential vanishes rendering 
the analytic derivation of quasi-normal frequencies impossible. We show that including the 
first-order correction leads to an analytic expression in agreement with numerical results. 
Unlike other types of perturbation, the correction in electromagnetic modes is logarithmic. 
We summarize our conclusions in section 5. 



2 Massive scalar perturbations 

In this section we calculate quasi-normal frequencies for massive scalar perturbations of 
finite black holes in AdS generalizing a procedure introduced in [2]. We consider explicitly 
the five-dimensional case in which the wave equation reduces to a Heun equation. General- 
izing to higher dimensions is straightforward albeit tedious due to the increase in singular 
points. 

Using the line element (1) in d = 5, we obtain the horizon radius 



r% 1 fl 2^7 

# = "2 + V4 + ^ (2) 



The wave equation for a massive scalar of mass m is 



f(r ) 



-^d r (r 3 f(r)d r $) - -7730*1 $ + — V 2 ^ = m 2 $ . (3) 



It is convenient to transform to a dimensionless coordinate 



V = s[ w + lj , * = *£^' W 
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in terms of which the factor f(r) (eq. (1) with d = 5) reads 

We see that s is a parameter describing the size of the black hole. When s — > 0, we approach 
the large black hole limit (rjj — ► oo) and expect to arrive at the results of [2]. 
Separating variables, 

$ = e-^Y m (n)^(y) , (6) 
we obtain the radial wave equation expressed in terms of y, 

~cj 2 (y - s) 2 L 2 



(y - s)(y 2 - + (3y 2 - 1 - 2sy) + 



4 y 2 - 1 4 
where we introduced the dimensionless parameters 



— - (y - s)m 



tf = (7) 



ijj 2 = 2su 2 R 2 , L 2 = 2s£(£ + 2) , m = — . (8) 

The singularities of the wave equation are given by 

y = ±l,s, (9) 

where y = 1 is the horizon, y = s is the black hole singularity and y = — 1 is an unphysical 
singularity. In order to bring (7) into a manageable form, we need to study the behavior of 
the wavefunction near the singularities. Two independent solutions of (7) are obtained by 
examining the behavior near the horizon (y — > 1), 

$r ± ~ (y - 1)±<|V^S . (10) 

where ^ r +,^ r _ represent outgoing and ingoing waves, respectively. We will choose \P_ for 
quasi-normal modes. 

Near the singularity y — > — 1 we obtain a different set of independent solutions 

^-(y+l) ± ^ vTTi . (11) 

Since this is an unphysical singularity there is no physical choice. By studying the behavior 
at large r (y — > oo), we find another set of independent solutions which determine the 
scaling behavior and are given by 



* ~ y~ h± , h± = l± ^/l + m 2 . (12) 

For quasi-normal modes we want the solution to vanish for large r (y — > oo), leading us to 
choose 

* ~ y~ h + . (13) 
We may write the solution of (7) in the form 

* = (y - + l)-I^F(y) . (14) 
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Substituting this expression into the wave equation (7), we obtain an equation for F(y), 
(y 2 - \)F" + |(3- {VTT~s + iVT- 



y + s + (VI + s - iVl-s)- 



F' 



+ 



s + iy 1 — s' 



{Vl + s + Wl- s) 



m 2 J> F 



+ - 



y- s 



(s 2 - l)F' -?±-F+[(l- s)VTT^ - i(l + s)v / T^] jF 



= (15) 

If we are interested in the limit of large frequencies u, we may focus on the region of large 
y [2]. In this case, the last term on the left-hand side of (15) is negligible compared with 
the other terms and the wave equation simplifies to a hypergeometric equation, 

(y 2 - 1)F" + 1(3- + y + s + (VTT~s - i^T^S)^-) f' 



s + i\J\- s 2 J - - (Vl + s + Wl - s) 



m 2 S> F 



(16) 



Two linearly independent solutions of (16) are 

T\ = F(a + , a_; c; —x) , Ti = x 1 ~ c F(l + a + — c, 1 + a_ — c; 2 — c; —x) 
where 

a± = h± — (Vl + s + iVl — s) 
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3 1/ • — ^ 

- + -(s-iVl-su) . 



(17) 

(18) 
(19) 



Using the transformation properties of hypergeometric functions, we may re-express the 
solutions (17) in terms of a new set of independent solutions which match the scaling 
behavior (12) for large r (x — > 00), 



K± = {x + l)- a± F(a±,c - a T ; a± - a T + 1; l/(x + 1)) 



(20) 



We ought to choose /C+, since it leads to \£ — > as x — > 00. K.+ may be expressed as a 
linear combination of T\ and J~2, 



where 



An 



T(l - c)r(l - a_ +a+) 



Bn 



r(c - i)r(i + a+ - a. 



r(l-a_)r(l-c + a+) r(a+)r(c-a_) 
For the correct behavior at the horizon, we demand 



(21) 
(22) 



B = 0, 



which leads to two conditions 



c — a_ = 1 — n , n = 1,2,3,... 



(23) 
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or 



a + = 1 — n , n = 1, 2, 3, 
Eq. (23) leads to the zeroth-order frequencies, 



-2(^l + s + Wl-s) 



n + h< 



(24) 



(25) 



Notice that the phase approaches 7r/4 in the large black- hole limit (r# — > oo or s — > 0), as 
expected [2]. 

Using (24), we find a second set of frequencies given by 



2(Vl + s - Wl-s)(n + h+ - 1) . 



(26) 



Both sets of frequencies, (25) and (26), at leading order agree on the imaginary part and 
have opposite real parts. We shall work with (25) without loss of generality. Notice also 
that the two sets of quasi-normal frequencies match the results of [2] in the large black hole 
limit (s -> 0). 

To find the first-order correction to the zeroth-order expression for quasi-normal fre- 
quencies (25), we shall solve the Heun equation (15) perturbatively. To this end, let us 
bring it to the form 

(H + H 1 )F = 0, (27) 

where (c/. eq. (16)) 







1 



y 1 y 2_ 1 



3 - (VTTs + n/T 



y + (s + {VTTs - iVl - s 



+- 



1 



U) 



(s + iy/l- s 2 )j - (Vl +s + Wl - s) 



V ~ 1 12 

and the correction (to be treated as a perturbation) is given by 
1 



m 



Hi 



(y 2 - i)(y - s) 



[s 2 - i)d y + ((i - s)VTT~s - i(i + s)vT^) 



On 



(28) 



(29) 



We have neglected the angular momentum contribution for simplicity. We may expand the 
wave function as 

F = F + Fi + . . . (30) 
where Fq obeys the zeroth-order equation (eq. (16)) 



n F = o. 



(31) 



Solving this equation leads to the zeroth-order expressions for quasi- normal frequencies (25). 
The first-order equation is 

TdFo + H F 1 = (32) 
We may solve for F\ by using variation of parameters, 



F, =/C_ 



K, + TL\Fq 

w 



IC^HiFq 

w 



(33) 



where K,± are the two linearly independent solutions (20) of eq. (16) and W is their Wron- 
skian given by 

W = (a+ -a_)x- c (l + x) c - a+ - a —\ (34) 
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To study the behavior near the horizon (x —►()), we may analytically continue the param- 
eters in (33) without affecting the singularity. For x ~ 0, we obtain 



F 1 ~A 1 + x 1 ~ c B 1 , 



where 



and 



fa 



JC_HiF 



With our choice (23), we find 



r(c-l)r(l + q ± -a T ) 
T(a±)T(c-a T ) 

K + HiF 



(35) 
(36) 
(37) 



Therefore the quasi-normal frequencies, to first order, are found as solutions of 

B + 0i = 0, (38) 

where Bo is given by (22). 

We can now find explicit expressions for the first-order correction to quasi-normal fre- 
quencies. Writing to first order 



LJ n = -2(Vl + s + iy/T 



3 s 
n + h + - 2 + 2~ £r - 



(39) 



we aim at calculating e n . Let us start with the case of n = 1. Our quantization condition 
(23) becomes c = a_ . This truncates the expansion of the hyper geometric solution (20) to 



After some algebra, we find 



Si 



fa 



2(a + — a_ 



F = /C + = (l + x)- a +. 



2>/„ dx ^TT^S 



k=0 



(40) 



(41) 



where the coefficients, a& (k = 0, 1), are given by 

a = -a+(s 2 -l) , ai= (s 2 - 1) +isy/l - s 2 [a+ - a_ + 1 + s] . (42) 

Using 

r°° t- a h 4- 

(43) 



o 



we find 
01 



00 x *(i _i_ x )-v 

dx^- =B(X + l,fi- X)F(1, A + 1; n + 1; 1 - 8), 

1 + ox 

B(a^ — 1, a + — a_ + 1 



1 - s 
a i 



a_|_ — a_ 



F(l, a_ + 1 a+ + 1; - 

s — 1 



2a 



s + 1 

— F(l,a_ + l,a+; -) 

2(a + — a_ — 1) s — 1 



(44) 
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Expanding in l//t+ (large mass expansion), we obtain 



Bi = B(a- - l,a+ - a_ + 1) 



| (s 2 - 1 + «sa/T 



8/l. 



(l + o( fl )) 



where we made use of the expansion of a hypergeometric function 



a 



+ 



1 1\ aV 



which is valid for large a and (3 (~ o(/i+)). We obtain from (22) and (39) 

Bo = e\B{a_ — 1, a + — a_ + 1) + . . . 
By using (38) we find the first-order correction for n = 1, 



+ 



"4 15 



1 + isVl-s 2 ) + ^- (1 + o(s)) 



8h-> 



(45) 



(46) 



(47) 



(48) 



For a finite-size black hole (s ^ 0), this is a c^/i 1 ^) correction to n = 1 quasi-normal 
frequencies. The correction is o(l//i_|_) for an infinite-size black hole (s = 0) [2]. It should 
be pointed out that the calculation of ei involved cancellation of o(h+) terms. For a general 
n, one obtains expressions o(/i"). Non-trivial cancellations occur between various terms 
involving hypergeometric functions and after the dust settles, one arrives at the general 
expression 



s 

An 



1 + is 



2(2 — l/n) , . . , , . 

+ g fe 7 ^ (l + o(a)) , 71=1,2,3,... (49) 



which is o(/i9_) for finite-size black holes and o(l//i+) for infinite-size black holes, 



1 



1 



1 



2n h 



(50) 



We have been unable to provide an analytical proof of the above results for general n but 
have verified them for several n using Mathematics. 



3 Gravitational perturbations 

In this section we discuss gravitational perturbations. For massless perturbations, the 
method discussed in section 2 is not directly applicable. Instead, we extend the proce- 
dure of [3, 4] to include first-order corrections to analytical expressions for quasi-normal 
frequencies. Our results are in good agreement with numerical results [7]. 

The radial wave equation for gravitational perturbations in the black-hole background (1) 
can be cast into a Schrodinger-like form, 



drl 

in terms of the tortoise coordinate defined by 

dr* 1 



+ F[r(r*)]* = cu 2 ^ , (51) 



dr f(r) ■ (52) 
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The potential V is determined by the type of perturbation and may be deduced from the 
Master Equation derived in [9]. For tensor, vector and scalar perturbations, we obtain, 
respectively, [4] 



£(£ + d -3) (d-2)(d-4)/(r) (d - 2)/'(r) 
Vr{r) = fir) { ^ 1 — h 



4r 2 



2r 



l(£ + d-3) (d - 2)(d - 4)/(r) rf"'(r) 
W(r-) = f(r) { 5 h 



4 r 2 



2(d - 3) 



(53) 
(54) 



Vs(r-) 



M 

4r 2 



+ d - 3) - (d - 2) + 



(d- l)(d-2)n 



„d—3 



d(d - l) 2 (d - 2)V 6(d - l)(d - 2) 2 (d - 4)[^(£ + d- 3) - (d - 2)]/i 



R 2 r 2d ~ 



R2 r d-h 



| (d - 4)(d - 6) [1(1 + d - 3) - (d - 2)] 2 r 2 + 2(d- l) 2 (d - 2)V 



i? 2 



r 



3d-9 



+ 



4(d - l)(d - 2)(2d 2 - lid + 18)[£(£ + d - 3) - (d - 2)]/x 2 



r 



2d-6 



+ 



(d - l) 2 (d - 2) 2 (d - 4)(d - 6)/i 2 6(d - 2)(d - 6)[£(i + d - 3) - (d - 2)] 2 ^ 



r 



2d-6 



6(d- l)(d-2) 2 (d-4)[^(£ + d-3) - (d - 2)]n 



„d~3 



+4[£{£ + d - 3) - (d - 2)] 3 + d(d - 2)[£{£ + d - 3) - (d - 2)] s 



(55) 



Evidently, the potential always vanishes at the horizon (V(r#) = 0, since f(ru) = 0) 
regardless of the type of perturbation. 

Near the black hole singularity (r ~ 0), the tortoise coordinate (52) may be expanded 

as 

^ yd— 2 i r 2d—5 

r * = ~{d-2) 2fi ~ (2d - 5) (2//) 2 + " ' (56) 

where we have kept the second term in the expansion of r and have chosen the integration 
constant so that r* = at r = 0. Using (56), we may expand the potential near the black 
hole singularity in the three different cases (eqs. (53), (54) and (55)), respectively as 



1 

4r 2 



+ 



A, 



-2(d-2)/x]3^3 



r* d " 2 + . 



At 



(d-3) 2 £(£ + d-3) 



2(2d - 5) 



+ 



d- 2 



Ay 



d — 1 
ti-2 



+ ... 



r * [-2(d-2)/i]3=J 



d 2 - 8d + 13 + d - 3) 
-4v = n , n , — — — I- 



2(2d- 15) 



d-2 



and 



where 



1 



-4s 



4r * [-2(d-2)fi]—2 



d-l 
d-2 



+ 



(2d 3 - 24d 2 + 94d- 116) (d 2 -7d+14)[^ + d-3)- (d-2)] 
-4s = 77^-, — — 1- 



4(2d-5)(d-2) 



(d- l)(d-2) 5 



(57) 
(58) 

(59) 
(60) 



We have included only the terms which contribute to the order we are interested in. We 
may summarize the behavior of the potential near the origin by 



V 



i' 2 — 1 ri-l 

+ Ar, d ~ 2 + ... 



4 r 2 



(61) 



where j = (2) for scalar and tensor (vector) perturbations and the constant coefficient 
A can be found from eqs.(57), (58), (59) and (60) in the various cases. Throughout the 
calculation, we shall pretend that j is not an integer. At the end of the calculation, we 
shall let j — ► 0, 2, as appropriate. 

After rescaling the tortoise coordinate [z = tor*), the Schrodinger-like wave equa- 
tion (51) with the potential (61) becomes 



d 2 ^ 
"cb 2 



+ 



Az 2 



1 



d-3 d-1 

-A oj d ~ 2 z d ~ 2 ^ 



(62) 



In the large frequency limit, we may treat the right-hand side of (62) as a correction. This 
will allow us to to solve the equation perturbatively. We may re-express (62) as 



where 



T~to = 



d 2 
dz 2 



Hq+UJ d-2 u x I * = 



4z 2 



, Hi 



d-1 

-Az rf- 2 . 



By treating Ti,\ as a perturbation, we may expand the wave function 

d-3 

V(z) = ^> {z) +u d - 2 ^\(z) + ... 
and solve (63) perturbatively. The zeroth-order wave equation, 

Ho^o(z) = 0, 
may be solved in terms of Bessel functions, 

^o(z) =A l ^fz~J ] (z)+A 2 y r zN i (z). 

2 2 

For large z, it behaves as 

^o(z) ~ [Aicos(z - a+) + A 2 sm(z - a + )] , 

= (A 1 - iA 2 )e- ia+ e lz + -J=(Ai + iA 2 )e +ia +e 
2vr 



(63) 
(64) 

(65) 
(66) 
(67) 



(68) 



where a± = \ (1 ± j). 

Next, we study the behavior of the wavefunction at large r. In this region, the tortoise 
coordinate (52) may be expanded as 



R 2 li? 4 

+ Q— + 

r Jr 



(69) 
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The integration constant is readily deduced from the definition (52) of the tortoise coordi- 
nate, 

(70) 



f(r) 

The potential (eqs. (53), (54) and (55)) for large r may be expanded as 

7 2 -1 

V= ., J °° _,, +••• (71) 

where joo = d — 1, d — 3 and d — 5 for tensor, vector and scalar perturbations, respectively. 
The Schrddinger-like wave equation (51) in the region of large r becomes 



+ 



n 2 - 1 



4(r* — 



* = (72) 



Since the potential does not vanish as r — ► oo, the wavefunction ought to vanish there. 
Imposing this boundary condition yields the acceptable solution to eq. (72), 



^(r*) = Byju)(r* - r*) J^u^r* - r*)) . (73) 

2 

Notice that <5 -> clS 7% ^ 7** , clS desired. Asymptotically, it behaves as 

*(r*)~l/f 5cos[w(r,-f*) + /3] , /3=~(l+joo) (74) 

By matching this expression to the asymptotic behavior (68) of the solution in the vicinity 
of the black-hole singularity along the Stokes line ^sz = ^(cjr*) = 0, we find a constraint 
on the coefficients A\, A 2 , 

Ai tan(wf* - (3 - a + ) - A 2 = 0. (75) 

A second constraint is obtained by imposing the boundary condition 

V(z) ~ e iz , z -> -00 , (76) 

at the horizon. To this end, we need to analytically continue the wavefunction near the 
origin to negative values of z. A rotation of z by —tt corresponds to a rotation by —3^2 
near the origin in the complex r-plane, on account of (56). Since near the origin, J u (z) ~ z v 
(multiplied by an even holographic function of z) and using the identity 

N v {z) = COtltU J v (z) — CSCTTV J- V (z) , (77) 

we deduce 

J v {eT iv z) = e-™ v J v (z) , N v {e~™z) = e™ v N v - 2i cos ttu J u (z) (78) 
Thus for z < 0, the wavefunction (67) changes to 

y (z) = e --(i+ 1 )/ 2 ^ | [Ax - i(l + e^)A 2 ] Ji(-z) + A 2 e™i N±(-z)} . (79) 
whose asymptotic behavior is given by 

* ~ ^— [At - i(l + 2e^ l ),4 2 ] e" 12 + =— [A x - iA 2 ] e lz (80) 
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Imposing the boundary condition (76) at the horizon, we deduce the constraint 



Ai - i(l + 2e j7Ti )A 2 = 0. 

The two constraints (75) and (81) are compatible provided 

1 + 2e j7ri ) 

tan(cjf* — P — a + ) —1 

which yields the quasi-normal frequencies [4] 

uf, = ~(2 + j + ioo) - tan" 1 - - ^ + nvr . 



0. 



(81) 



(82) 



(83) 



These are zeroth-order expressions deduced from the zeroth-order wave equation (66). 

Next, we calculate the first-order correction to the asymptotic expressions (83) for quasi- 
normal frequencies. We begin by focusing on the region near the black-hole singularity 
(r ~ 0). To first order, the wave equation (63) becomes 



where Hq and 7i\ are given in eq. (64). The solution is 

^ Ji(z')HiMz') 



(84) 



9 1 (z) = VtN i (z) / dz' 



>V 



z Jj (z) / dz' 



z> NjWHiVoW 

_2 

w 



(85) 

written in terms of the two linearly independent solutions (67) of the zeroth-order eq. (66). 
W = 2/-7T is their Wronskian. Using (67) and (85), we may express the solution to the wave 
equation (63) up to first order (eq. (65)) explicitly as 



= {Aiil - b(z)\ - A 2 a 2 (z)} J~zJ 3 (z) + {A 2 [l + b(z)] + Aioi(z)} J~zN, {i 

2 2 

where the functions a\(z), a 2 (z) and b(z) are given by 

■kA 



ai(z) 



d-3 
CO ' d-2 



, , irA 
a 2 {z) = —^-u d ~ 2 



b(z) 



IT A d-3 

2 



dz' Z f -7=3J i (z')J i (z f ), 
2 2 

dz' z'~^N L (z')N L (z), 

2 2 



(86) 



(87) 



(89) 



respectively. The coefficient A is defined in eq. (64) and depends on the type of perturbation. 
The wavefunction (86) behaves asymptotically as 



[2 , , 
ty(z) ~ d- [A 1 cos(z — a + ) + A 2 sm{z — a+)] , 



where 



A[ = [1 - b]Ai - a 2 A 2 , A' 2 = [1 + b]A 2 + a\A x 
and we introduced the notation 

d\ = oi(oo) , a 2 = a 2 (oo) , b = b(po) . 



(90) 
(91) 
(92) 
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By matching this to the asymptotic expression (74), we obtain 

A 1 tan(cjf* — (3 — a+) — A 2 = 



(93) 



correcting the zeroth-order constraint (75). Using (91), the first-order constraint (93) in 
terms of A\ and A 2 reads 

[(1 - b) tan(wf» -/?-«+)- ai]At - [1 + b + a 2 tan(u;f» - (3 - a+)]A 2 = (94) 

To find the first-order correction to the second constraint (81), we need to approach the 
horizon. This entails a rotation by — ir in the z-plane. Prom the small- 2 behavior of a Bessel 
function, J u (z) ~ z" ', and using the identity (77), we deduce after some algebra 



ai (e- in z) 
b(e~ in z) 



e -^ e -^ai(z) 
e d ~ 2 



2 



e m3 a 2 (z) - 4cos" yai(z) - 2i(l + e^)6(z) 



e -^— [6( z ) _ j(i + e -^) ai (^) 



(95) 



From these expressions and eq. (78), we arrive at a modified expression for the wavefunc- 
tion (86) valid for z < 0. In the limit z — > -co, we obtain 



-le 



-wPBx cos(-z - a+) - ie ijn / 2 B 2 sm(-z - a+) 



(96) 



where 



B l = A l -A 1 e~™— [b - i{l + e-^cn] 



-A 2 e 



e +inj a 2 - 4 cos 2 H gl - 2i(l + e +i ^)b 



5 2 = A 2 + A 2 e- in ^[b-i(l + e- i7Tj )a 1 ]+A 1 e- iw ^e- i7Tj ai 
By imposing the boundary condition (76) at the horizon, we obtain 



(97) 



[1 _ e -'*a=5 (ia! + 6)] Ai - [i(l + 2e^') + e" m ^ ((1 + e i7rJ >i + e^a 2 - ib)]A 2 = (98) 

correcting the zeroth-order constraint (81). For compatibility of the two first-order con- 
straints, (94) and (98), we need 



1 + b + a 2 tan(wf* - (3 - a+) i(l + 2e™i) + e~ i7r ^((l + ^)a x + e^a 2 - i&) 
(1 — b) tan(a;r* — (3 — a+) — a\ 1 — e d ~ 2 [iax + b) 

Solving (99), we arrive at the first-order expression for quasi- normal frequencies, 



= 
(99) 



7T 1 

Lor„ = -(2 + j + joo) + — In 2 + nil 
4 2% 



If- _,vM=3 T -iir d - 3 -iir d - 3 

--\6ib-2ie J7r ^6-9ai + e m ^a 1 +a 2 - e m d-2 
8 L 



«2 



(100) 
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where we took the limit of interest j — ► 0,2 wherever it was unambiguous, in order to 
simplify the notation. Using 



r(A)r( ^+ 1 ~ A ; 



d,X X X J u (x)J u (x) = ; — —r rTT r T~ ! rTH , (101) 

A*V 1 "\ J \ -m-i / — iy-t-//-t-1 4-i\xi / ?y — //-t-1 -t-\\-n/iy-l-;(-t-14-X \ ' V / 



-i^+/l+l + A -jp^ t/-/J+l + A ^-p ^+^t+l + A 

we obtain explicit expressions for the first-order coefficients, 

ttA / nir\~i^ r (d=2") r (2 + 2(d-2)) 

ai 



2 



2r *' ^ ( 2{d-2) 2 + 2(d-2) ■ 

, „ 7r(d- 3) vr / . d- 3 
1 + 2 cot -f- f cot — [ —J + 



2{d-2) 2 V d-2 



01 



I = -cot ^ 3 j at (102) 
2(d-2) v ; 

where we used the identity r(rr)r(l — x) = si ^ nx ■ We also set to = mr/f*, since corrections 
contribute to higher than first order. Notice that these expressions are well-defined when j 

d-3 

becomes an integer. Thus, the first-order correction is ~ o(n d ~ 2 ). 

Next, we compare with numerical results in four dimensions [7]. It is convenient to set 
the AdS radius R = 1. From (1), the radius of the horizon r# is related to the black hole 
parameter fi by 

2fi = r 3 H + r H (103) 
for d = 4. f(r) has two more (complex) roots, r_ and its complex conjugate, where 

r = e'W3 + - e — (J^2 + JL + ^ (104 ) 

The integration constant in the tortoise coordinate (70) is 

00 dr r_ , r_ r* r* 

111 — -^2— T ln — (105) 







jf(r) 3ri +1 r# 3r* 2 + 1 r# 



Despite appearances, this is not a real number, because we ought to define arguments as 
< argr < 2ir. 

For scalar perturbations, we find from eqs. (100), (102), together with (59), (60) and 
(61), 



Notice that only the first-order correction is ^-dependent. In the limit of large horizon 
radius (th ~ (2/i) 1 ^ 3 S> 1), we have from (105) 



_ _ 7r(l + Z>/3) 
3V3r H 

Numerically for I = 2, 



(107) 



uj n , , 0.508 + 0.293* 
— = (1.299 - 2.250i)n + 0.573 - 0.419i H 108) 

r H rfj ' 
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Figure 1: The frequency gap (110) for scalar perturbations in d = 4 for r# = 1 and £ = 2: zeroth 
and first order analytical (eq. (109)) compared with numerical data [7]. 
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Figure 2: The frequency gap (110) for tensor perturbations in d — 4 for r h — 1 and £ = 0: zeroth 
and first order analytical (eq. (114)) compared with numerical data [7]. 



For an intermediate black hole, rjj = 1, we obtain 

uj n = (1.969 - 2.350i)n + 0.752 - 0.370i + °- 654 ± ( 1Q9 ) 



In figure 1 we compare this analytical result with numerical results [7]. We plot the gap 

Aw n = w n -w n _i (110) 

because the offset does not always agree with numerical results [4]. We show both zeroth- 
order and first-order analytical results. For a small black hole, rjj = 0.2, we obtain 

1.093 + 0.561t , x 

oj n = 1.695 - 0.571»)n + 0.487 - 0.044H + = (111) 

In 



For tensor perturbations, we find from eqs. (100), (102), together with (57) and (61), 
^. , S-o , ^/4^TF 4 (i) rj^ 3£(£ + l) + l 



" + i)-^ 2 -w-ii_^_, At ^^^i (112) 
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Figure 3: The frequency gap (110) for tensor perturbations in d = 4 for r# = 0.2 and I 
zeroth and first order analytical (eq. (115)) compared with numerical data [7]. 
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Again, only the first-order correction is ^-dependent. Numerically for large m and £ = 0, 

^ = (1.299 - 2.250*> + 0.573 - 0.419, + °- 102 + ;_ 0586 * (113) 
r H rfjy/n 

For an intermediate black hole, rjj = 1, we obtain 

u n = (1.969 - 2.350i)n + 0.752 - 0.370i + °- 131 + 0-0916? 

\ n 



In figure 2, we plot the gap (110), including both zeroth and first order and compare with 
numerical results [7]. 

For a small black hole, rjj = 0.2, we obtain 

u n = (1.695 - 0.571t> + 0.487 - 0.044H + °' 489 ± J' 251 * ( 115 ) 



and compare the gap with numerical results in figure 3. Finally, for vector perturbations, 
we find from eqs. (100), (102), together with (58) and (61), 



n+ 4r + 2 ln2+e ^8^v^' - 4v= ^ — + n (116) 



Numerically for large th and I = 2, 

^ = (1.299 - 2.250i)n + 0.573 - 0.419* + 8A9 + 6 ^ dl ( 117 ) 

For an intermediate black hole, rjj = 1, we obtain (see figure 4) 

uj n = (1.969 - 2.350i)n + 0.752 - 0.370i + °- 741 + 0-519? 

and for a small black hole, rn = 0.2, we obtain (see figure 5) 

1 239 ~\~ (33572 

uo n = (1.695 - 0.571*)n + 0.487 - 0.0441* + = (119) 

\/n 



In all cases of gravitational perturbations, regardless of the size of the black hole, our 
analytical results are in good agreement with numerical results [7]. 
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Figure 4: The frequency gap (110) for vector perturbations in d = 4 for r# = 1 and £ = 2: zeroth 
and first order analytical (eq. (118)) compared with numerical data [7]. 
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Figure 5: The frequency gap (110) for vector perturbations in d — 4 for r# = 0.2 and 
zeroth and first order analytical (eq. (119)) compared with numerical data [7]. 
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4 Electromagnetic perturbations 

In this section we extend the discussion to electromagnetic perturbations. This is a singular 
case because the potential vanishes at zeroth order. Consequently, the compatibility con- 
dition (82) discussed in the previous section has no solutions and no asymptotic expression 
for quasi-normal frequencies may be deduced [4]. Nevertheless, the numerical results are 
similar to the ones we discussed in the case of gravitational perturbations [7], We shall 
show that including first-order corrections leads to analytical asymptotic expressions for 
quasi- normal frequencies in agreement with numerical results. Unlike with gravitational 
perturbations, where first-order corrections were a power of n (eqs. (100) and (102)), for 
electromagnetic perturbations first-order corrections are o(lnn). 

We shall concentrate on the four-dimensional case for defmiteness. Generalization to 
higher dimensions is straightforward. The wave equation reduces to (51) with electromag- 
netic potential 

Vem = ^J^/(r). (120) 
where f(r) is given in (1) with d = 4. Near the origin, this potential may be expanded in 



16 



(122) 



terms of the tortoise coordinate. Using eq. (56), we obtain 

FEM "^T + 2^ + --" (121) 

where j = 1. This leads to a vanishing potential to zeroth order. Consequently, no analytic 
expression for quasi-normal frequencies is deduced. This is easily seen by substituting j ' = 1 
in the zeroth-order expression (83); we obtain a divergent result because tan -1 i is not finite. 

This is remedied by including first-order corrections. The compatibility condition (99) 
of the two first-order constraints (94) and (98) reads 

1 + b + a,2 tan wf* —i — b + iai 
(1 — b) tan wf* — Oi 1 — a\ + ib 

where we used d = 4, j = 1, j'qo = d — 3 = 1, and a + = (3 = \ . At zeroth order (setting 
o-i = 0,2 = b = 0), we obtain 

tanuf* = i (123) 
which has no solution, as expected [4]. At first order, we obtain 

tanwf* = i + (1 -i)(dx -a 2 - 26) (124) 

whose first-order solution is 

1 (l + i)(ai-q 2 -25) 

o;r* = rnr H In (125) 

2i 2 

Using (102) and (121), we deduce explicit expressions for the first-order coefficients, 

ai = V^> a 2 = 6 = -ai, -4=^4= (126) 
and eq. (125) reads explicitly 

i 1 

cjf* =nvr- -Inn + — In (2(1 + ^^v 7 ^) (127) 
4 2i 

Therefore, the correction to the quasi-normal frequencies behaves as Inn in the large lo 
limit. 

To compare with numerical results, set R = 1. As with gravitational perturbations, we 
shall compare the gap, because the offset is not reliable. For the gap, we have from (127) 

AWn = UJ n - UJ n -l = y (l - ^ + . . )j (128) 

Both leading and sub-leading terms are independent of t. 

For a large black hole, using (107), we obtain the spectrum 

Aw„ 3^/3(1-^^/3)/ * \ 0.179 + 0.103i . s 

m y — Z—L i _ — + . . . = 1.299 - 2.25i + . . . 129) 

th 4 y Ann ) n 

This analytical result is compared with numerical results [7] for th = 100 in figure 6. 

Using eqs. (103), (104), (105) and (127), we obtain the spectrum of an intermediate 
black hole, m = 1, (see figure 7) 

u n = (1.969 - 2.350i)n - (0.187 + 0.1567i) In n + . . . (130) 

and for a small black hole, th = 0.2, (see figure 8) 

uo n = (1.695 - 0.571z)n - (0.045 + 0.135t) In n + . . . (131) 

All first-order analytical results are in good agreement with numerical results [7]. 
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Figure 6: The frequency gap (110) for electromagnetic perturbations in d = 4 for r# : 
6=1: zeroth and first order analytical (eq. (129)) compared with numerical data [7]. 
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Figure 7: The frequency gap (110) for electromagnetic perturbations in d = 4 for m = 1 and 
£ = 1: zeroth and first order analytical (eq. (130)) compared with numerical data [7]. 




Figure 8: The frequency gap (110) for electromagnetic perturbations in d = 4 for r# = 0.2 and 
£ = 1: zeroth and first order analytical (eq. (131)) compared with numerical data [7]. 
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5 Conclusions 



We studied quasi-normal modes for Schwarzschild black holes in asymptotically AdS spaces 
of arbitrary dimension. We obtained analytical expressions by solving the wave equation 
perturbatively, including first-order corrections. We studied scalar, electromagnetic and 
gravitational perturbations. In the case of massive scalar perturbations, we extended the 
method proposed in [2] for large black holes and derived explicit expressions of quasi-normal 
frequencies for black holes of arbitrary size as a perturbative expansion in 1/m, where m 
is the mass of the perturbation. 

This method is not directly applicable to massless modes, because the perturbative 
expansion fails as m — > 0. Instead, we obtained perturbative expansions of gravitational 
and electromagnetic modes by extending the method proposed in [6] for asymptotically flat 
spaces. The perturbative expansion was based on zeroth-order results obtained in [3, 4]. We 
showed that our analytical results were in good agreement with numerical data [7]. In the 
case of electromagnetic perturbations, zeroth-order expressions do not yield finite quasi- 
normal frequencies, because the effective potential vanishes [4]. By including first-order 
effects, we were able to arrive at finite analytical expressions with logarithmic sub-leading 
contributions. 

It would be interesting to extend these results to other types of black holes and also 
understand their implications on the AdS/CFT correspondence. 
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